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ABSTRACT 

This paper presents new observations of the AGNs M87 and Hydra A at 90 GHz made with 
the MUSTANG bolometer array on the Green Bank Telescope at 8.5" resolution. A spectral 
analysis is performed combining this new data and archival VLA data on these objects at longer 
wavelengths. This analysis can detect variations in spectral index and curvature expected from 
energy losses in the radiating particles. M87 shows only weak evidence for steepening of the 
spectrum along the jet suggesting either re-acceleration of the relativistic particles in the jet or 
insufficient losses to affect the spectrum at 90 GHz. The jets in Hydra A show strong steepening 
as they move from the nucleus suggesting unbalanced losses of the higher energy relativistic 
particles. The difference between these two sources may be accounted for by the different lengths 
over which the jets are observable, 2 kpc for M87 and 45 kpc for Hydra A. 

Subject headings: galaxies: jets, galaxies: active, radio continuum, galaxies: individual (M87, Hydra A), 
GBT 



1. Introduction 

It has long been recognized that the radiative 
lifetimes of the synchrotron emitting electrons in 
AGN jets and extended radio lobes can be shorter 
than the time to transport them from the nu- 



cleus to wher e 
IWillis fc Strg 



they 



are observed ( Felten 19681 : 
19781 ). This requires some resup- 



ply or re-acceleration of the radiating particles 
well outside of the nucleus. Observ ations of syn- 



chrot ron jets visible into the X-ray (jBiretta et al 



^also University of Maryland, College Park, MD 20742- 
2421 



199ll ) make this problem particularly acute. For 
a synchrotron source of a given magnetic field 
strength, the emission observed at a given fre- 
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quen cy comes largely fr om electrons of a giyen en- 
ergy (|Pacholczyklll970l ). Thus, the spectral shape 
of the emission reflects the shape of the electron 
energy spectrum. If the predominant energy loss 
mechanism of these electrons is radiative (syn- 
chrotron and inverse Compton scattering) then 
the losses are proportional to the energy squared. 
This depletes the number of electrons above a cor- 
responding energy resulting in a steepening of the 
spectrum by 0.5 above a break frequency corre- 
sponding to the highest energy surviving particles 
(^acholczyk 1970,). If the radiating electrons age 
in a constant magnetic field, the break frequency 
will move to progressively lower frequencies. 

Numerous attempts have been made to use this 
spectral signature to determine the radiative ages 
of the particles, i.e. the time since t hey were last 



accel e rated to relativistic energies (ICarilli et al 



1991 



Blundell &: Alexandeil 11994 iFeretti et al 



1998 1: 'Murgia et a.1.1 119991 ). As was pointed out 



by iRudni ck (2003) , this procedure has a number 
of potential difficulties. Most serious is that the 
relationship of the frequency of the spectral break 
with the cutoff in the electron energy spectrum 
depends on a large po wer of the poorly known 
magnetic field strength; Pacholczva ( 1970?) gives 
the relation among the break frequency(;/B), the 
magnetic field(i?) and time(t) as: 



const B t 



1) 



Thus a lower magnetic field results in higher break 
frequencies. While using spectra to determine par- 
ticles ages is suspect, the steepening of the spec- 
trum along an AGN jet, especially at higher fre- 
quencies shows the relative effects of particle losses 
and gains. 

Most of the sources in which the signatures of 
spectral aging have searched for have been ob- 
served over a relatively limited frequency range. 
Since the effects of the electron energy losses are 
more pronounced at higher frequencies, observa- 
tions at high frequencies are valuable for evalu- 
ating the relative effects of particle losses against 
any re-acceleration or resupply. 

In order to measure the variation of the ob- 
served spectrum along an AGN jet, both good 
resolution and surface brightness sensitivity over 
a wide range of frequencies is needed. Observa- 
tions with mm interferometers have the resolution 
but generally not the surface brightness sensitiv- 



ity for comparisons with lower frequency observa- 
tions. Since single dishes explicitly sample all spa- 
tial frequencies out to the dish diameter, they have 
good surface brightness sensitivity. The large size 
of the Green Bank Telescope (GET) gives it good 
resolution as well as surface brightness sensitivity. 

In this paper we report on the observations of 
M87 and Hydra A at 90 GHz using the MUS- 
TANG camera on 100 m t he Green Bank Telescope 
( Jewell fc Prestagj 2004 ) which gives 8.5" resolu- 
tion and good surface brightness sensitivity. These 
results are compared with lower frequency obser- 
vations at comparable resolutions using archival 
VLA data. The results presented give the con- 
tinuum spectra from 0.3 to 90 GHz; this greatly 
increases the range of frequencies over which this 
analysis can be performed. 

2. Observations 

The radio sources in the AGNs M87 and Hy- 
dra A were observed on 2008 March 13 using the 
MUSTANG 3mm bolometer camera (jPicker et al 



20081 ) on the Green Bank Telescope (GBT). The 
MUSTANG detector array consists of 64 TES 
bolometers viewing the sky through a filter with 
bandpass of 81 to 99 GHz. Observations were 
made in the "On-The-Fly" (OTF) mode using a 
variety of scanning patterns. A cold field stop re- 
stricts the portion of the GBT illuminated to 90 m. 
The bolometer outputs are recorded at 1 kHz but 
the data streams are averaged to 20 Hz for anal- 
ysis. Focus and pointing calibration used J1256- 
0547 for M87 and J0825+0309 for Hydra A. Flux 
density calibration was based on observations of 
planets made over the course of the season's ob- 
servations. 

3. MUSTANG Data Analysis 

Dat a analysis wa s carried out in the Obit pack- 
age (jCottonl 2008f l. Details of the data analysis 
will be given in a future publication but a short 
description follows. 

The principle difficulty in imaging the data 
from the MUSTANG array is separating the por- 
tion of the signal in the time-stream due to the 
astronomical object of interest and the "back- 
ground" portion due to other sources, especially 



-^http: / /www. cv. nrao.edu/~bcotton/Obit. html 
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the atmosphere and the instrument itself. Rapid 
scanning of the telescope over the target field helps 
separate the signals in the time domain. The re- 
dundancy resulting from the multiple observations 
of a given pointing along different trajectories also 
helps separate the time variable background sig- 
nals from the constant astronomical signals. 

The MUSTANG cryogenic system uses a pulse 
tube cooler with a frequency of 1.4 Hz which im- 
poses a strong modulation on the bolometer out- 
puts. Prior to subsequent operations, the ampli- 
tude and phase of the resultant signal were fitted 
and removed from the data. Amplitude calibra- 
tion used the response to a stable calibration sig- 
nal injected into all bolometers whose strength was 
calibrated against observations of several plan- 
ets. Data weights were derived for each bolome- 
ter time-stream from the RMS fluctuations during 
blank sky scans. 

3.1. Imaging 

Imaging used the general "OTF" technique 
of Mangum. Emerson, fc GreisenI (|2007h in which 
the individual detector samples were multiplied 
by a continuous "griding" function and then re- 
sampled on a regular grid with 2" spacings. For 
each grid cell, the weighted sum of the data times 
the griding function and the sum of the weights 
times the griding function were accumulated and 
the image value at that pixel was the set equal to 
the ratio of the two. The griding function used 
was an exponential x sine which helps suppress 
spatial frequencies larger than those sampled by 
the dish and results in a resolution of 8.5". This 
process will produce a "Dirty Image", i.e., an im- 
age of the sky convolved with the telescope beam 
response. 

3.2. Background Estimation 

An iterative scheme of estimating the target 
brightness distribution, subtracting a model from 
the time-stream data and (re)estimating the back- 
grounds was employed. In each iteration the short- 
est timescale used to estimate backgrounds was 
reduced. 

For the initial estimation of the backgrounds 
an empty sky was assumed and backgrounds were 
estimated for timescales longer than 40 seconds. 
This is longer than the crossing times of the tar- 



get sources. Following the initial background es- 
timation, and for subsequent iterations, the back- 
ground subtracted time-stream data were imaged 
as outlined in Section [37T] A deconvolved image 
(see § 13. 4p was used to generate a sky model which 
was then used with the assumed telescope point 
spread function (PSF) to generate a residual time- 
stream data set from the initial calibrated time- 
streams. A refined background estimation was ob- 
tained from a low pass filtering of the time-stream 
residuals. This refined background estimation was 
in turn used to derive a refined sky model. Multi- 
ple iterations of this process were used, gradually 
reducing the shortest time scale in the background 
estimations to 1 second. Most of the backgrounds 
are common to all detectors. 

3.3. GBT Imaging Quality at 90 GHz 

At the epoch of these observations, uncorrected 
deformations of the GBT surface caused a degra- 
dation of the telescope performance. In particu- 
lar, large scale errors in the surface result in an 
elevated error beam in the neighborhood of the 
main beam. While the near-in error beam is time 
variable in detail, general features persisted over 
time. 

If the point spread function (PSF) of the instru- 
ment is know and stable, it is relatively straight- 
forward to deconvolve the dirty image and replace 
it with one in which the psf was well behaved. The 
best estimate of the effective PSF of the telescope 
was obtained from the averaged image of calibra- 
tor observations made during the 2008 Winter ob- 
serving session. A "Dirty beam" for the telescope 
was derived by an azimuthal average of the aver- 
age calibrator beam images below the half power 
of the beam and an 8.5" Gaussian above the half 
power point. This beam has an elevated error 
pattern near the main beam and when used in a 
deconvolution, removes the scattered power from 
around regions of bright emission. Since these 
data were collected, a method of applying online 
corrections to the active surface derived from out- 
of-focus beammaps has significantly reduced the 
scattered-power problem. 

3.4. Deconvolution 

The dirty image and dirty beam described in 
the previous sections were deconvolved using a 
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Hogbom CLEAN (jHogboml 120031 ) . This reduces 
the dirty image to a set of delta functions which 
can be convolved by the GBT instrumental PSF 
to produce a sky model or convolved with a Gaus- 
sian and restored to the CLEAN residuals to get 
the final CLEAN image. These CLEAN images 
restored with a 8.5" FWHM Gaussian were used 
in subsequent analysis. 

One consequence of the way in which the dirty 
images were made is that noise in the image pixels 
does not strictly have the same covariance as the 
dirty beam. This allows "out of band" noise which 
does not deconvolve and causes the CLEAN to 
"get stuck", alternately digging and filling holes. 
A filtering of the dirty image to remove all power 
corresponding to spatial frequencies on smaller 
scales than allowed by the resolution of the in- 
strument greatly reduces this problem. 

The images of M87 and Hydra A derived in the 
Obit package are given in Figures [T] and 

4. VLA Data at Longer Wavelengths 

In order to make a proper spectral analysis, im- 
ages at comparable resolution are needed over a 
range of frequencies. Fortunately, these bright ra- 
dio galaxies have been well studied using the VLA. 
Suitable resolution data were obtained from the 
VLA archives and calibrated and imaged in the 
standard fashion using Obit. The imaging was ini- 
tially at 8" resolution using a combination of uv 
ran ges, tapering , and the Briggs robust parame- 
ter ( Briggd 1995[ ). Subsequently, the images were 
convolved to 8.5" for comparison with the GBT 
data. 

5. Spectral Decomposition 

The GBT and VLA images were all aligned 
to the same pixel grid and two sets of spectral 
forms were fitted in each pixel by standard least 
squares techniques. These forms were a contin- 
iously curved spectrum with polynomials in log(s) 
vs log(i^) using up to three terms; 1) the flux den- 
sity at a reference frequency, 2) the spectral index 
and 3) spectral curvature at the same reference 
frequency. The other form was a broken power 
law where the spectral index at frequencies above 
the break steepened by 0.5. All fitting used a ref- 
erence frequency of 5 GHz. For each source, a set 
of single point spectra were derived as well as im- 
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Fig. 1.— Cleaned Obit MUSTANG image of M87 
at 8.5" resolution. Contours are powers of ^/2 
times 0.1 Jy/beam. 
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Fig. 2.— Cleaned Obit MUSTANG image of Hy- 
dra A at 8.5" resolution. Contours are powers of 
\f2 times 0.025 Jy/beam. 
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ages of the spectral index, spectral curvature and 
break frequency. 

Plots of the derived image and spectral param- 
eters are shown in Figures [3] for M87 and [5] for 
Hydra A. Letters on Figures [3^ and [5^ give the lo- 
cations where individual spectra were determined. 
These individual point spectra are in FiguresHand 
[6] for M87 and Hydra A respectively. Both contin- 
uously curved spectral fits and broken power laws 
are shown. Corresponding fitted parameters are 
given in Tables [1] and [2l These tables give the 
flux density in Jy, spectral index and curvature 
terms (when fitted), psr degree of freedom of 
the curved spectrum fit, the break frequency in 
GHz from the broken power law fit and ^ for the 
broken power law fits at the locations indicated in 
Figures [4] and [6l 

6. Discussion 

The presence of optic al and X-ray emis sion 
from the inner jet in M87 (jBiretta et al.lll99ll ) ne- 
cessitates a resupply of energetic particles in at 
least this portion of the jet. The data presented 
above can probe the possibility of further particle 
resupply in the portions of the jet only detected 
in radio emission. 

The radiative "age" of a source with break fre- 



quency oivB is given bv lCarilli et al.l (jl99ll ) as 

t ^ f6fOS-i-^z/^°-^ 2) 



where B is in micro-gauss, vb in GHz and t in Myr. 
Thus, the electrons radiating in a f micro-gauss 
field giving rise to 90 GHz radiation will start to 
be seriously depleted in about 5 Myr. 

6.1. M87 

One of the earliest cases recognized to require 
particle acceleratio n outside of t he nucleus was the 
optical jet in M87 (|Feltenlll968h . M87 is the dom- 
inant elliptical galaxy in the center of the nearby 
Virgo cluster. Previously p ublishe d radio obser- 
vation s inc lude Owen et al. (f980), iBiretta et al.l 
(|f983l ) and lOwen et all (|2000( ). We assume a dis- 
tance to M87 of f 6 Mpc which results in a linear 
projected size of 78 pc per arc-second. The jets are 
thus visible out to beyond 2 kpc from the nucleus 
in Figure [3l 

M87 is an asymmetric two-sided FRI jet source 
in which both jets bend into the line of sight lim- 
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Fig. 3.— a) Fitted image of M87 at 5 GHz with 
indications of locations of point spectra as letters. 
Gontours are powers of \/2 times 25 mJy/beam. 
b) Grey-scale spectral index with contours of flux 
density at fOO, 300, 1000 and 3000 mJy/beam 
over-plotted. The gray-scale is given by the scale- 
bar at the top. The resolution is indicated by the 
circle in the lower left corner, c) as b) but break 
frequency (GHz) shown in gray-scale. 
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Fig. 4. — Point spectra of M87 at locations indicated in Figure [H Fitted values are shown in Table [T] Solid 
lines show the curved spectrum fits and the dashed lines broken power law fits. 
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Fig. 6. — Point spectra of Hydra A at locations indicated in Figure [5l Fitted values are shown in Table [2l 
Solid lines show the curved spectrum fits and the dashed lines broken power law fits. 
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Table 1: M87 Fitted Point Spectra at 5 GHz. 



Posn. 




Flux 


Spectral 


Curvature 


A curve 


Break Freq. 


-^break 




density (Jy) 




index 










GHz 




A 


3.9 ± 0.2 


-0.62 ± 0.03 


-0.05 ± 0.01 


1.5 


3.3 ± 1.2 


2.2 


B 


6.1 


0.3 


-0.56 


0.03 


-0.05 


0.01 


1.8 


3.0 


0.2 


2.2 


C 


4.9 


0.2 


-0.29 


0.03 


0.03 


0.01 


1.4 


3.9 


0.8 


13.2 


D 


2.6 


0.1 


-0.74 


0.02 


-0.04 


0.01 


0.3 


3.9 


0.1 


3.3 


E 


4.1 


0.2 


-0.67 


0.03 


-0.04 


0.01 


0.7 


3.4 


0.1 


3.5 


F 


2.7 


0.1 


-0.65 


0.03 


-0.05 


0.01 


1.3 


3.3 


0.1 


2.0 



Table 2: Hydra A Fitted Point Spectra at 5 GHz. 



Posn. 




Flux 


Spectral 


Curvature 


A curve 


Break Freq. 


'^break 




density (Jy) 




index 










GHz 




A 


0.4 ± 0.03 


-1.20 ± 0.08 


-0.11 ± 0.03 


4.7 


6.4 ± 0.8 


5.4 


B 


3.2 


0.22 


-0.80 


0.03 


-0.07 


0.01 


0.1 


4.9 


0.1 


0.2 


C 


0.4 


0.05 


-0.47 


0.07 






1.1 


8.2 


3.2 


0.4 


D 


1.6 


0.11 


-0.71 


0.03 


-0.04 


0.01 


0.6 


4.1 


0.1 


2.7 


E 


0.8 


0.05 


-0.83 


0.03 






3.0 


8.3 


4.4 


1.04 



iting the distance from the nucleus to which they 
can be observed. Figure [3] contains the jets and 
inner lobes embedd ed in the mu c h larg er region of 
emission shown bv I Owen et al.l ( 2000l ). The jets 
on either side of the nucleus (position C) show lit- 
tle variation in spectral index and curvature. The 
approaching jet (positions A, B) have somewhat 
flatter spectra than in the receding jet (positions, 
D, E, F). The fitted break frequencies shown in 
Figures [3] and [4] and Table [1] show minimal vari- 
ation along the visible portion of the jets. Fur- 
thermore, the continiously curved spectral model 
gives consistently better fits, in terms of the 
per degree of freedom than the broken power law 
spectra. 

Magnetic field estimates based on minimum en- 
ergy arg uments and higher resolution VLA images 
given bv lHines et al. ( 1989f ) for features in the re- 
gions probed here, range from 50 to 100 micro- 
gauss. This leads to estimates of radiative life- 
times of the 90 GHz emitting electrons in the range 
150 to 500 kyr. We note that this is easily within 
the possible transport times for the radiating par- 
ticles from the nucleus. 

The variation in spectral index and curvature 
across the source shown in Figures |3] and 2] and 
Table [T] show very weak steepening with distance 



from the nucleus. The curved nature of the spec- 
tra do indicate high energy losses to an assumed 
initial power law electron energy spectrum. The 
better fit of the curved spectrum than a simple 
broken power law model suggests a complex parti- 
cle loss/resupply history rather than a simple, ra- 
diatively evolving relativistic particle population. 
Since the observable portion of the jet is relatively 
short, the lack of evidence for electron aging could 
be due to either continued particle acceleration 
along the jet or merely insufficient time for the 
particle populations to evolve to the point they 
could be detectable in our data. 

The asymmetry of the spectral index is likely 
due simply to the emission on the counterjet side 
being stongly contaminated by the steeper spec- 
trum "cocoon". The jet, being Doppler boosted, 
dominates the cocoon on its side. 

6.2. Hydra A 

The radio structur e of Hydra A has been stud- 
i ed b y I Taylor et al.l 
(Il99,5[ ) and iLane et al 



1990|) .Dwarakanath et al 



2004[ l: extended regions 



of very steep spectrum emission are seen to the 
north and south of the emission in Figure O This 
source is a symmetric two-sided FRI radio galaxy 
with a weak, strongly self absorbed core (po- 
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Fig. 5. — a) Fitted image of Hydra A at 5 
GHz with indications of locations of point spec- 
tra as letters. Contours are powers of ^/2 times 
10 mjy/beam. b) Grey-scale spectral index with 
contours of flux density at 10, 30, 100 and 300 
mJy/beam over-plotted. The gray-scale is given 
by the scale-bar at the top. The resolution is indi- 
cated by the circle in the lower left corner, c) as b) 
but break frequency (GHz) shown in gray-scale. 



sition C). Hydra A is at a redshift of z=0.054 
( Dwarakanath et al. |[l99i) giving a projected lin- 
ear size of 1.1 kpc per arc-second; the source is 
visible in Figures [H and [S] to approximately 45 kpc 
from the nucleus. 

Magnetic field estimates of lTavlor et al. I (|l990l) 
from minimum pressure arguments range from 
about 10 to 70 micro-gauss in various features giv- 
ing estimated radiative lifetimes of the 90 GHz 
emitting electrons in the range 0.3 to 5 Myr. 

In this galaxy, both the spectral index and 
curvature increase dramatically along the visi- 
ble portion of the jets, especially in the north. 
This suggests that particle acceleration slows or 
stops fairly early in the jet and that the effects 
of particle aging are visible in the outer parts 
of the sourc e . Th e magnetic field estimates by 
Taylor et al. ( 1990[ ) suggest weaker field strengths 
in the northern lobe. By Eq. 1, weaker fields 
sho uld reduce the effect of electro n aging (however, 
see iBlundell fc Alexandeil [l99l . In the north- 
ern part of the source (locations A and B) the 
curved and broken spectra are comparable repre- 
sentations of the observed spectra whereas in the 
south (location E), the broken power law spec- 
trum more nearly represents the high frequency 
measurements. However, in spite of strong varia- 
tions in spectral index, there does not appear to 
be a strong variation in the break frequency as 
represented in Figures [S] and [S] and Tabled if any- 
thing, the break frequency increases slightly with 
distance from the nucleus. This is more sugges- 
tive of a decreasing magnetic field than an aging 
electron population. The apparent break in the 
spectrum in the south may be the result of syn- 
chrotron losses in the inner regions of the jet where 
the magnetic field is expected to be stronger. 
This is inconsistent with a simple intrepretation 
of the steepening along the jet being due to elec- 
tron losses; the fitted break frequency may not 
be meaningful except possibly in the south where 
the broken power law is a better representation of 
the spectrum than the curved model. The curved 
spectrum in the north suggests a complex electron 
history. 

The observable portion of the jet in Hydra A is 
substantially longer than in M87. The longer time 
required to transport the radiating particles to the 
outer regions allows radiative losses to accumulate. 
This may be responsible for the apparent break in 



the spectrum in the south. 

7. Conclusions 

We present new 90 GHz observations of the FRI 
AGN radio sources M87 and Hydra A at 8.5" res- 
olution and compare the results with lower fre- 
quency archival VLA results at comparable reso- 
lution. A spectral analysis is performed to look for 
evidence of particle aging or re-acceleration. M87 
shows only weak variation out to 2 kpc but Hydra 
shows strong steepening and bending over regions 
extending to 45 kpc. This difference may be the 
result of the different size scales observed in the 
two sources. In M87, the length of jet observed 
may be too short for either the particle resupply 
to damp out or for the radiating particles to age to 
the point the the effects are visible in the observed 
spectral region. A continiously curved spectrum 
gives consistently better fits than a broken power 
law spectrum. This is more indicative of a complex 
history of the electron energetics than a simple ag- 
ing of a population of radiating particles. 

In Hydra A there is clear evidence of a deple- 
tion of higher energy electrons, the loss mecha- 
nisms are winning over the resupply mechanisms. 
Even in the case of Hydra A, there is little evi- 
dence for a variation in the break frequency with 
distance from the nucleus. Only in the southern- 
most observed portion on the jet is there credible 
evidence of a break in the spectrum. The gener- 
ally curved nature of the observed spectra suggest 
complex histories of particle loss and resupply. 
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